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Abstract. By means of analytical calculations, we explore the self-enrichment scenario for Globular Cluster forma- 
tion. According to this scenario, an initial burst of star formation occurs inside the core radius of the initial gaseous 
distribution. The outward-propagating shock wave sweeps up a shell in which gravitational instabilities may arise, 
leading to the formation of a second, metal-enriched, population of stars. We find a minimum mass of the proto- 
globular cluster of the order of 10 6 Mq . We also find that the observed spread in the Magnitude-Metallicity relation 
can be explained assuming cluster-to-cluster variations of some parameters like the thermalization efficiency, the 
mixing efficiency and the Initial Mass Function, as well as variations of the external pressure. 
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1. Introduction 

Globular Clusters (GCs) are among the oldest objects in 
the Galaxy and thus are fossil records of its formation. 
Their age is comparable with the age of halo stars but, 
quite surprisingly, the lowest metallicity found in GCs 
([Fe/H] < —2.5) is more than two orders of magnitude 
larger than the iron abundance of the most metal-poor 
halo stars (e.g. the recently discovered star HE 0107-5240 
has a metallicity [Fe/H] = -5.3; Christlieb et al. 2004). 
This could be only a statistical effect, since GCs make up 
only ~ 1% of the stars in the halo of the Galaxy. Moreover, 
the metallicity distribution of the halo stars and of the 
GCs peak at the same value (Ashman & Zepf 1998). This 
could mean that the halo stars are the debris of the tidal 
disruption of cluster-like objects. Observations of the dis- 
tribution of a wide sample of stars detected in the Sloan 
Digital Ski Survey (Newberg et al. 2002) seems to support 
this picture. 

More detailed analyses of the metallicity distribution 
of halo stars and GCs (Laird et al. 1993; Carney et al. 
1996) indicate that there is a statistically significant dif- 
ference between the two metallicity distributions that can- 
not be explained by the lower numbers of GCs. Moreover, 
field halo stars are almost all CN-weak, whereas some GCs 
show a bimodal CN distribution, with some stars having 
high CN abundances (Norris & Smith 1981). This may 
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indicate a difference in the formation process of these two 
classes of objects. 

In many GCs the abundances of elements from C to Al 
show complex patterns and star-to-star abundance varia- 
tions, such as the well-known O-Na and Mg-Al anticor- 
relations (Gratton et al. 2001; Gonzalez & Wallerstein 
1998). A few GCs: uj Centauri (Freeman & Rodgers 1975; 
Smith et al. 2000; Pancino et al. 2002), M22 (Norris & 
Freeman 1983; Lehnert, Bell & Cohen 1991) and maybe 
M92 (Langer et al. 1998) show a significant spread even 
in iron-peak abundances, with signatures of a metallicity 
gradient. Pancino et al. (2002) found also in u> Centauri 
that [Fe/H] correlates with [Cu/Fe], whereas the [Si/Fe] 
vs. [Fe/H] plot shows a plateau until [Fe/H] ~ —1 and 
then an anticorrelation. These patterns are characteristic 
of iron enrichment from SNela (Matteucci et al. 1993). 

The range of metallicities in GCs is similar to that 
found in Local Group Dwarf Spheroidals (DSphs) and 
even in this case no galaxies with [Fe/H] < —2.2 are found. 
This may indicate a common origin of these two classes of 
objects, or at least a similar enrichment scenario. 

In order to explain the lack of metal-poor GCs, two 
processes have been invoked: pre-enrichment from PopIII 
stars (Beasley et al. 2003) and the self- enrichment sce- 
nario, the second being the only one able to explain abun- 
dance variations observed in some GCs. The idea behind 
this is that, if a proto-globular cluster (PGC) is massive 
enough, it may retain the heavy elements produced by 
the first supernovae exploding after the collapse of the 
cloud. These metals may mix with the surrounding, unpol- 
luted medium and any following episodes of star formation 
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arise from a metal-enriched gas. The pioneers in this kind 
of study have been Dopita & Smith (1986) and Cayrel 
(1986), but many other authors have studied the subject 
(Morgan & Lake 1989; Ikuta & Arimoto 2000; Smith 2000 
among others). In particular, Shustov & Wiebe (2000; 
hereafter SW00) and Parmentier et al. (1999; hereafter 
P99) addressed the question of how massive a PGC should 
be in order to retain the products of at least one supernova 
(SN) and how many SNe can a PGC sustain before being 
disrupted as a result of the injected energy. P99 tried also 
to compare the results of their model with the magnitude- 
metallicity relation in a particular sample of Milky Way 
GCs (namely the older ones). 

This is the first of a series of three papers dealing with 
analytical and numerical calculation of the self-enrichment 
process in GCs. In this paper, after carefully studying 
the available data about GCs and Local Group DSphs 
(Section 2), we perform an analytical study of the evolu- 
tion of a PGC, taking into consideration the possibility 
of triggered star formation (Section 3). Results are pre- 
sented and discussed in Section 4. Finally, we draw our 
conclusions in Section 5. In the forthcoming papers, we 
will analyze 1-D (paper II) and 2-D (paper III) chemody- 
namical simulations of the evolution of a PGC. 
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Fig. 1. Plot of the total visual magnitude My vs. metal- 
licity [Fe/H] for GCs of the Milky Way (squares) and 
Local Group DSphs (triangles). Data are from the Harris 
Catalogue (GCs) and from Mateo (1998) (DSphs). 



2. On the mass dependence in the metallicity 
distribution of Milky Way GCs 



The most obvious reference 
of Milky Way GCs is the 



for any study 
Harris Catalogue 



I http://physun.physics.mcmaster.ca/Globular.html 
see also Harris 1996). The metallicity-luminosity diagram 
for the whole galactic globular cluster system is shown in 
Fig. 1. We include also in the same plot the metallicity- 
luminosity relation for Local Group DSphs, taken from 
the Mateo (1998) review. 

This plot does not show any obvious correlation among 
GCs, whereas the DSph system exhibits a well defined and 
well known trend, such that the faintest DSphs are also 
the most metal-poor. The classical explanation of this cor- 
relation is that the effect of galactic winds decreases with 
the increasing depth of the galactic potential well (Larson 
1974; Dekel & Silk 1986). This mechanism cannot easily be 
extended to GCs, mostly because the DSphs show a very 
large mass-to-light ratio, presumably a hint of dominant 
dark matter halos, whereas GCs show very little evidence 
of dark matter. 

There have been several attempts in the literature to 
calibrate the integrated light of GCs in order to recover 
the metallicity (Danziger 1973; Zinn 1980). The Harris 
catalogue is based on the Zinn & West (1984) calibration, 
but this is not the only one and care is needed before 
drawing conclusions from the inspection of the metallicity- 
luminosity plot. In particular, Carretta & Gratton (1997) 
proposed a revised version of the Zinn & West (1984) 
calibration. This new calibration changes the resulting 
metallicity in particular in the range of high abundances. 
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Fig. 2. Plot of the total visual magnitude My vs. metal- 
licity [Fe/H] for GCs of the Milky Way, by means of the 
calibration proposed by Carretta & Gratton (1997). 



The metallicity-luminosity relation with this calibration is 
shown in Fig. 2. 

It has been suggested that the mass loss due to strip- 
ping with passage through the disk of the Galaxy might 
dominate the mass budget of GCs. This could be the rea- 
son why there is such an enormous spread in this relation. 
One may suppose that the galaxies with large thus 
at large distance from the disk of the Galaxy, may have 
suffered less stripping. For this reason, we select the GCs 
with 1 6 1 > 20 and plot in Fig. 3 the metallicity-luminosity 
relation for this selected sample of GCs, with both the 
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Fig. 3. Plot of the total visual magnitude My vs. metal- 
licity [Fe/H] for the GCs of the Milky Way with |6| > 20. 
The filled squares are relative to the Zinn & West (1984) 
calibration, the open squares to the Carretta & Gratton 
(1997) one. 



Zinn & West (1984) calibration (filled squares) and the 
Carretta & Gratton (1997) one (open squares). 

Although Fig. 3 shows a reduced spread compared to 
Figs. 1 and 2, an obvious correlation still cannot be de- 
tected, suggesting that, besides external effects, internal 
properties of the GCs should also be invoked in order to 
justify this enormous spread. 

For a reduced sample of GCs, mass measurements are 
available. In this case, we do not have the problem to 
calibrate the metallicity index, but the determination of 
the mass of GCs is very uncertain and model-dependent 
(Pryor & Meylan 1993). We adopt the same selection cri- 
terion introduced for Fig. 3 (i.e. selection of GCs at large 
galactic latitudes) and we plot in Fig. 4 the GC masses as 
a function of the metallicity for the GCs studied by Pryor 
& Meylan (1993). Also in this case, no obvious trend can 
be identified, unless we exclude from the plot the two most 
massive GCs. In this case we can see that the metallicity 
slightly increases with mass. 

3. An analytic approach 

3.1. Set-up of the model 

Let us assume an isothermal sphere (P ~ p) in hydrostatic 
equilibrium. The hydrostatic equilibrium, in spherical 
symmetry, is expressed by the formula ^ = — GM ^ p ^ ; 
whereas the conservation of mass implies that =p- = 
4irp(r)r 2 . Combining these two equations, we obtain the 
fundamental equation of equilibrium: 



l_d_ 

r 2 dr 



r 

p dr 



(1) 



Fig. 4. Plot of the total GC mass log Mqc vs. metallicity 
[Fe/H] for the GCs of the Milky Way with \b\ > 20. Data 
are taken from Pryor & Meylan (1993). 



We can modify this equation, introducing two dimension- 
less variables 



u = - In p/po 

e = Ar, A 2 = ^P-, 



where c s is the sound speed, obtaining: 



2 dui 



(2) 



(3) 



This is the classical Lane-Emden equation for an isother- 
mal sphere and cannot be solved analytically. With a 
simple, second-order numerical integration, we obtain the 
density profile plotted in Fig. 5. Also plotted (dotted line) 
is the asymptotic trend at large £ (p ~ po ■ £~ 2 ) 

This density profile can be well approximated by a 
King profile: 



p{r) 



Po 
Po 1 



if r < r c 
if r > r c , 



(4) 



where r c is the core radius and the profile is truncated 
at a tidal radius R. The concentration parameter c = 
log(i?/r c ) is ~ 1 for most of metal-poor GCs (see Harris 
Catalogue), thus we can consider hereinafter r c — 0.1 R. 
We assume pressure equilibrium between the PGC and 
the external medium, thus, defining Ph as the external 
pressure, we obtain p(R) = 0.01 po — Phjc 2 . We use this 
expression to calculate the mass of the PGC as a function 
of R and Ph . Finally, this expression can be used to obtain 
the tidal radius R as a function of the total PGC mass and 
the external pressure. The final expression is: 



R 



1 3/cT 
2.8 AnfirriH 



1/3 




1/3 



(5) 
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Fig. 5. Normalized density profile in an isothermal 
sphere. Solid line: solution of the numerical integration; 
dotted line: asymptotic trend at large £. 

Assuming /j, = 1.3 and T = 10 4 (as we will do hereinafter), 
we obtain: 



R 100 = 1.385 f ^ 



1/3 



(6) 



where i?ioo is in units of 100 pc, Mq in units of 10 6 M Q 
and P4 in units of 10 4 fc. The dependence of po with Ph is 
Po,22 = 2.I7P4, where po,22 is in units of 10 -22 g cm -3 . 

3.2. The collapse of the PGC 

Let us now assume that the sphere loses its equilibrium 
and starts collapsing. This assumption is justified by the 
fact that in the eq. (1) we neglected the cooling, thus 
the configuration should lose the equilibrium due to ra- 
diative losses. The cooling time in the inner core (t c ~ 
3kT / (2nA(T)) , where A(T) is the cooling function) is ap- 
proximativcly 2 -10 4 yr with our assumptions. The collapse 
of the core will be much faster than the collapse of the 
external mantel. If we assume for simplicity the free fall 
time-scale of the mantel as tff. man — ^/[37r/(32Gp man )], 
where p man = 3po?" 2 {R — r c ) / (R 3 — r 3 .) is the average den- 
sity of the external mantel, we obtain £//, c A//,mtm ~ 0.16; 
that is, if we assume for simplicity two separate collapse 
episodes, the mantel will not have time to loose its equi- 
librium configuration in the interval of time in which the 
core collapses. We can thus assume that a fraction of the 
core collapses and forms stars in an instantaneous burst of 
star formation, while the external mantel keeps its starting 
density profile. 

The efficiency of star formation in the core (namely 
the fraction of core mass turned into stars) is assumed 



to be 0.3, consistent with the star formation efficiencies 
obtained in the cores of dense clouds by Bate, Bonnell & 
Bromm (2003). 

The mass of stars in this first, almost metal-free, 
episode of star formation is thus: 



1.071 • 10" 2 M fi 



(7) 

The number of SNell is given by: Nsnii = f$ ° Am~ x dm, 
where A is calculated from M* = f 100 Am 1 ~ x dm. A 
Salpctcr (x = 2.35) index can be taken as a good ap- 
proximation to the IMF slope in clusters (Kroupa 2002). 
mi is the lower cut-off of the IMF. We define g(mi) = 
100" 35 and we take mi as a free parameter. We 



m 



-0.35 



obtain: 



N s 



Nil 



162.1- 



M fi 



(8) 



9 {mi) 

For a reasonable choice of mi (i.e. mj = 1) this gives 
Nsnii ^ 200 M 6 , thought by P99 to be the criti- 
cal number of SNe sustainable against PGC disruption. 
Assuming the rate of SNell constant in time, we can 
calculate the luminosity of the burst. This is given by: 
L = tjEsnii Nsnii/ At, where r\ is the fraction of the me- 
chanical energy of the explosion which goes into thermal 
energy of the medium (the so-called thermalization effi- 
ciency), Esnii is the energy of a single SNII (assumed to 
be 10 51 erg) and At ~ 30 Myr is the interval of time in 
which SNell explode. The thermalization efficiency can be 
evaluated assuming the evolution of a SNR in a uniform 
medium. We can calculate the time at which the expan- 
sion velocity of the SNR becomes equal to the local sound 
speed, the kinetic energy of the shell at this time, and 
divide this energy by the blast wave energy (see Cioffi, 
McKee & Bertschinger 1988; Bradamante, Matteucci & 
D'Ercole 1998; Recchi, Matteucci & D'Ercole 2001 for 
more details). We obtain: 



n 



0.02 no 54/49 r 15/9% cJl 



(9) 

where £ is the metallicity (in units of Z©) and co,6 is 
the local sound speed (in units of 10 6 cm s _1 ). Assuming 
£ = 10~ 4 (typical of the most metal-poor halo stars), we 
obtain: 



r] = 4.36 • 10" 4 P 4 



-54/49 



(10) 



It is worth mentioning that the determination of 77 in this 
way is too simplistic and many other physical processes 
contribute to its determination (see Melioli & de Gouveia 
Dal Pino 2004). Therefore we will keep r\ in this paper as a 
free parameter, bearing in mind that, at least in the initial 
stages of the evolution of the system, its value might be 
very low. The luminosity of the burst (in units of 10 38 erg 
s _1 ) is thus given by: 



L 38 = 1.712M 6 



9{mi) 



(11) 



We assume that mass and energy of this first generation of 
stars are released in a spherical volume of radius tsf « 
R and we keep the r~ 2 gas density profile outside it. 
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Table 1. Ejecta fraction and mass in metals 
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2.34 
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2.74 
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0.93 


4.06 


30 


0.94 


6.01 


35 


0.92 


7.18 


40 


0.90 


8.22 



Total ejected fraction (f e j — 1 - M rem /Mi„) and mass of metals 
m z as a function of initial mass. Data from Woosley & Weaver 
1995 (Z=10~ 4 Z ; case B) 

3.3. Slow winds or fast winds? 

The evolution of wind-blown bubbles and supcrbubbles 
has been carefully analyzed by Koo & McKee (1992a; 
1992b). Since the results of these authors will be exten- 
sively used in the paper, we recall here briefly a few re- 
sults about bubble expansion. The freely expanding wind 
produced by the starburst interacts supersonically with 
the unperturbed ISM, creating a classical bubble struc- 
ture (Weaver et al. 1977), in which two shocks are present. 
The external one ("ambient shock") propagates through 
the ISM and creates an expanding cold and dense shell. 
A second shock ("wind shock") propagates inwards, ther- 
malizing the impinging wind and creating the hot, rarefied 
gas of the bubble interior. The shocked ambient medium 
and the shocked wind are separated by a contact disconti- 
nuity. Initially, the wind density is large and the velocity 
of the wind shock is small, therefore the wind shock is ra- 
diative. At later stages, the wind density decreases and the 
wind shock velocity increases, therefore the wind shock is 
expected to become adiabatic. If the energy input is pow- 
erful enough, this transition could occur in an early stage 
of the bubble expansion, when the wind is still freely ex- 
panding. In this case, the radiative cooling of the cavity 
is dynamically unimportant. Koo & McKee (1992a) call 
these kinds of winds fast winds. The slow winds are in- 
stead the ones in which the transition from a radiative to 
an adiabatic wind shock occurs later, making the cooling 
of the bubble dynamically relevant. 

Considering a generic power-law density distribution 
of the ambient medium p{r) = poir~ kp and a generic rate 
of energy injection L in = Ct^' 1 , Koo & McKee (1992b) 
were able to calculate a critical velocity, namely the veloc- 
ity which divides the two types of winds. This expression 
is given by: 



( L 



\ 27n 7i. 



3p, 



'01 



3 k r 



2-77; 



1 



ci 



3-fc„-7? irl 



,(12) 



where ctk, v = 1/[14 — 6fc p — (3 + fc p )?7i„] and C\ = 6 x 10 



-35 



g cm 6 s 4 . 



Assuming a r 2 profile and a constant injection rate 
(i.e. 7]i n = 1) cq. (12) becomes: 



C 



1/3 



6irp r 2 



(13) 



By substituting our parameters, we obtain: 

1/9 / \ 1/3 



v cr , 8 = 2.84 • 10" 




V 



(14) 



An approximate value of the velocity of the ejecta from 
SNell is given by 

E. 



= \/2E SNII /M ej , SN ii, where 
snii — V ■ Nsnii ■ 10 51 erg is the total energy released 
by SNell, whereas M e j^Nii is the total mass restored by 
these SNe. This mass is a fraction f e j of the total mass 
in stars in the interval [8, 100] M Q . Assuming a Salpeter 
IMF, the mass of stars in this interval is 0.283 / g{m{) 
times the total mass of stars, whereas f e j can be evaluated 
from the tables of Woosley & Weaver (1995). Assuming an 
abundance of the PGC of 10~ 4 Z Q , f e j <~ 0.9, with almost 
no dependence on the initial mass (see Table 1), we thus 
obtain: 



2.440?, 



(15) 



where v e j t g is in units of 10 8 cm s 1 . 
The ratio v e j/v cr is given by: 



85.9 




-1/9 



1/3 



(16) 



which is very large and varies very little with the param- 
eter space. Even assuming the (very low) thermalization 
efficiency found in cq. (10), assuming m; = 1 and neglect- 
ing the very weak dependence on M 6 and P4, we obtain 



20. We can thus safely state that the winds 



occurring in PGC are fast winds. 



3.4. The evolution of the Superbubble 

Koo & McKee (1992a) defined the fiducial radius Rf as 
the radius at which the wind density equals the ambient 
density. For a general power-law density distribution and 
a general rate of energy injection, they obtained: 



(3 - k p )C 



(17) 



Koo & McKee (1992) were able to derive the expansion 
law of the bubble under various initial conditions and in 
different stages of the bubble evolution, as a function of 
Rf. Indeed, ours is a special case, since, for k p — 2 and 
rji n = 1 this normalization breaks down and the fiducial 
radius approaches zero. This special case is called con- 
stant velocity bubble and, physically, this breakdown cor- 
responds to the fact that the ratio between swept-up mass 
and ejected mass is constant. 
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The ambient medium always dominates these kinds of 
winds and the bubble always stays radiative. The super- 
bubble expands with a constant velocity = v e jl/(l + l), 

where / = <J £ / (2ir po^v^) = V3 ■ (v cr /v ej ) 3 / 2 . With our 
parameters: 

1/6 



v b6 ~ 0.532 




(18) 



(where ti^g is in units of 10 6 cm s 1 ) and consequently 
the superbubble expands with a law: 



R 



1/6 



■s.lOO 



P 



(19) 



with t§ in units of 10 6 yr. Since the burst originates in a 
spherical region of radius r c , R s should be considered as 
the difference between the shock radius and r c . Note that, 
if the bubble evolves into a partially radiative bubble (i.e. 
a bubble in which most of the shocked wind has cooled 
down, but most of the bubble volume is filled with the 
most recently shocked and still hot portion of the wind), 
the velocity would be no more constant, but increasing 
with t 1 / 3 . However, such an accelerating bubble would 
produce Ray lcigh- Taylor instabilities with the overlying 
gas, resulting in a larger cooling of the shocked wind, pre- 
venting the onset of the partially radiative bubble (see Koo 
& McKee 1992b). Therefore, the assumption of a constant- 
velocity bubble is reliable. The time needed to reach the 
tidal radius R is: 



= 22.9 



Me 
Pi 



1/6 



V 



-1/4 



9 (mi) 



1/2 



(20) 



which is of the order of the lifetime of a 8 M Q star. It 
implies that the assumption of a constant luminosity is 
reliable. 

The swept-up mass as a function of time is given by: 



m sw , 6 = 4.2 • io- 2 M 6 5 ^ /f y /4 s(™ ; r 1/2 *6. 



(21) 



It is easy to check that the swept-up mass at tud is similar 
to the initial mass; the difference being the mass in stars. 

The mass mz of metals ejected in the ISM by a SN 
whose progenitor mass is Mi n is shown in Table 1. This 
dependence can be approximated by: 



m z (M in ) = 0.015(M m - 8) : 



H • 10- 5 (M in - 8) 4 , (22) 



(see Fig. 6). This approximation is valid until Mj„ = 40 
Mq. We will assume a constant mz above 40 Mq. The 
total mass of metals released into the ISM as a function 
of time can be expressed as: 



M z (t) 



0.35 

9 (mi) 
0.35 
9 (mi) 



A'L 



i(t) 



M- I - ib m z (M m )dM vt 



M*A z (t). 



(23) 



0.04 



0.02 




5 10 15 20 25 
t - t(100) (Myr) 



Fig. 6. Normalized metal production Az as a function of 
time. r(100) is the lifetime of a 100 M Q star. In the inner 
plot is shown the mass mz of heavy elements ejected by 
a SN of initial mass Mi n and our approximation to this 
trend. 



This integral cannot be evaluated analytically. We solve 
it numerically. We assume a simple law for the stellar life- 
times: t(m) = 1200m" 185 -3 Myr (Padovani & Matteucci 
1993). The resulting Az(t) is shown in Fig. 6 and the val- 
ues every Myr are tabulated in Table 2. 

3.5. Triggered star formation in the shell 

The expansion of a supershell can trigger star formation 
in the shell itself if the onset of gravitational instabili- 
ties in the fragments occurs when the shell is still bound 
to the PGC potential. Moreover, we must require that 
the time-scale for the onset of these gravitational insta- 
bilities is shorter than the time-scale needed to reach the 
tidal radius of the PGC. In order to calculate the time- 
scale needed for the onset of gravitational instabilities, 
we follow an approach similar to the one of McCray & 
Kafatos (1987), namely we consider a small circular disk 
with radius £ < < r on the surface of the expanding spher- 
ical shell, which expands in a r~ 2 density profile with 
a constant velocity U&. We indicate with 8 ~ £/r the 
half-angle subtended by the disk. The thermal, kinetic 
and gravitational energy of the disk are mc 2 . im,(-S-} 



s ' 2 



. dt I 



nQ 2 por 2 {r — r c ), thus we obtain: Ek 



, . . , , respectively. The mass of the disk is m(r,9) = 

±Tr6 4 v 2 p rl(r-r c ); 



E Th = ir6 2 p r z (r- 



„2. 



Eb 



GTT 2 8 3 p 2 4(r-r c ) 2 /r. 



The criterion for the onset of a gravitational instability is 
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Table 2. Normalized metal production Az as a function 4. Results 
of time. 

4.1. Constraints on the mass of the PGC 
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Normalized metal production as a function of time (see eq. 21 
for the definition of Az(t)). 



approximativcly Ek 
into: 



Etk + Eb < 0, which transforms 



+ c 2 < 0. 



(24) 



The A of this equation must be > 0, otherwise the above 
expression is positive for any value of 9. This transforms 
into the following condition: 

-l 

3c s v b 



r > r c 



1 - 



V2Gnp r 2 c 



(25) 



provided that 3c s Vb < \/2GTrp r 2 . 



t*,6 



This radius is reached at a time 
2.55(M f ) 1/ V 1 / 4 g(mi )i/2 



1 - 11.07(M 6 P 4 )- 1/2 ?7 1/4 9H) _1/2 ' 



(26) 



The metallicity of this self-enriched population of stars is 
simply given by: 



Z = 1 M Z (U)/M SW (U), 



(27) 



where 7 is the mixing efficiency, namely the fraction of 
metals produced in the first generation of stars, able to 
mix with the surrounding cold shell in a time-scale shorter 
than the time-scale needed for the onset of gravitational 
instabilities in the shell. 



The time needed for the onset of a gravitational instability 
able to generate the population of stars we observe nowa- 
days in GCs (t*^) is given by eq. (26). This time-scale has 
to be shorter than the timc-scalc needed to reach the tidal 
radius (t t id,6, given by eq. (20)), otherwise the energy re- 
leased by the first generation of stars simply blows away all 
the gas in the PGC. By comparing these two time-scales, 
we obtain a condition for the mass of the PGC: 



M M »n,i > 12.46F 4 ^n^girmY 



1/2 



(28) 



We can obtain another constraint on M 6 by considering 
that, in the moment in which stars form in the shell, the 
shell has to be bound. It means that the potential energy 
of the shell should be larger than its kinetic energy, i.e. 



lGM, 
2 " 



,(t,y 



minimum shell mass: 



In this way we can define a 



M 



shminM — 



9.H. lQ~-2^) 2 l\l/2 g{mi) -l 

1 - 11.07(M 6 P 4 )- 1/2 r? 1 / 4 3(m / )- 1 /2- 



(29) 



The swept-up mass calculated by means of eq. (21) has to 
be larger than M s h m in in order that the shell be bound. 
This translates into another condition for the mass of the 
PGC, namely: 



M Mm . 2 > Q.&P^ 2 rf /2 g{ mi )- 



(30) 



As usual, all the masses (eqs. 28, 29 and 30) arc in units 
of 10 6 M . 

The resulting dependence of the minimum mass of the 
PGC calculated according to eq. (28) (Mmm) and eq. 
(30) (MMin,2) is shown in Fig. 7 (as a function of n) and 
in Fig. 8 (as a function of mi). For simplicity, the external 
pressure is kept constant in these plots and in the following 
sections. It is assumed to be P4 = 10, in agreement with 
the pressure of the hot protogalactic background (Fall & 
Rees 1985; Murray & Lin 1992). We analyze the depen- 
dence of [Fe/H] on P4 in Sect. 4.2.4. These plots show 
that for any reasonable choice of parameters, the critical 
mass required in order to get the shell bound to the PGC 
potential well (M Mi „ i2 ) is always extremely low. This is 
therefore not a severe constraint on the minimum mass of 
the PGC. 

The constraint based on the time needed for the on- 
set of a gravitational instability compared to the time re- 
quired for the shock to reach the tidal radius (expressed 
by eq. (28)) is more stringent and led us to consider only 
PGC with masses larger than ~ 10 6 M . This is the 
same threshold mass found by SW00 and is also compa- 
rable with the threshold mass found by Brown, Burkcrt 
& Truran (1995), but is significantly larger than the one 
found by Morgan & Lake (1989). 
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log (?;) 




Fig. 7. Lower limits for the mass of the PGC calculated 
according to eq. (28) (upper panel) and eq. (30) (lower 
panel). The thermalization efficiency 77 varies with conti- 
nuity, whereas the lower IMF mass mi is equal to 0.1 (solid 
line), 1 (dotted line), 3 (dashed line). 

4.2. The Magnitude-Metallicity relation in GCs 

A direct comparison between the results of this simplified 
model and the observations can be made as follows. We 
can assume the same efficiency of star formation adopted 
in Sect. 3.2, namely that 30% of the mass of the shell 
at the moment of the onset of the gravitational instability 
(M sw (t*)) transforms into stars. This is therefore the mass 
of our newly formed GC. This should be considered as 
a lower limit, since, after the formation of the second, 
polluted generation of stars, the GC is still rich in gas. 
Shocks created by the supernovae exploding in the shell 
can still compress the gas, leading to the formation of new 
stars. 

To convert the mass in stars into visual magnitude, 
we assume a constant mass-to-light ratio M/Ly = 1.5, 
typical for GC systems in the Milky Way (Harris 1996), 
M31 (Djorgovski et al. 1997) and M33 (Larsen et al. 2002). 
At this point the visual magnitude is simply calculated by 
means of the following formula: 



5.41 — 2.5 • log [0.2 • M sw (t*)] . 



(31) 



Now we can directly compare the observed magnitude- 
metallicity relation with the results of our approxima- 
tions. We take the sample of data already shown in Fig. 3, 
namely the GCs of the Milky Way with \b\ > 20, in order 
to avoid the effect of the mass loss due to stripping passing 
through the disk of the Galaxy. We explore the parameter 
space, varying the mixing efficiency 7, the thermalization 



Fig. 8. As Fig. 6 but with varying lower IMF mass mi. 
The thermalization efficiency r\ is equal to 10~ 5 (solid 
line), 10~ 3 (dotted line), 0.1 (dashed line). 



efficiency 77, the lower limit of the IMF distribution mi 
and the external pressure P4. 

4.2.1. The Magnitude-Metallicity relation as a function 

of 7 

The mixing between the metals produced in the first 
generation of stars and the swept-up shell is a complex 
process, involving poorly known physics. In the idealized 
wind-blown bubble model, the shocked wind (hot, cen- 
tral region filled with metals) and the shocked ISM (the 
shell of unpolluted swept-up material) are separated by 
a contact discontinuity, thus in principle no flux of mat- 
ter is possible between these two regions. Actually many 
physical processes can allow a mixing (mixing layers, tur- 
bulence, thermal instabilities, thermal conduction), but 
these processes are poorly constrained. We therefore de- 
cided to keep the mixing efficiency 7 as a free parameter, 
trying to find some constraint on its value through the 
comparison of our results with the observations. 

The Magnitude-Metallicity relation as a function of 7 
is shown in Fig. 9. In this set of models, the thermaliza- 
tion efficiency 77 and the lower limit of the IMF distribu- 
tion mi are set to be 10 -4 and 1, respectively. Although 
the value of the thermalization efficiency is rather low, it 
is fully consistent with the value found in eq. (10). As we 
can see in this figure, the minimum required mixing effi- 
ciency is between 0.05 and 0.1, whereas, in order to reach 
the most metal-rich clusters, a very large 7 is required. A 
non-negligible mixing between the hot cavity, filled with 
metals, and the shell, should therefore occur, but the mag- 
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Fig. 9. Observed Magnitude-Metallicity relation for a se- 
lected sample of GCs (see Fig. 3) compared with model 
results. The mixing efficiency 7 is taken as 0.01 (solid line), 
0.05 (dotted line), 0.1 (dashed line), 0.5 (long-dashed line) 
and 1 (dot-dashed line). 



Fig. 10. Observed Mass-Metallicity relation for a selected 
sample of GCs (see Fig. 4) compared with model results. 
The mixing efficiency 7 is taken as 0.01 (solid line), 0.05 
(dotted line), 0.1 (dashed line), 0.5 (long-dashed line) and 
1 (dot-dashed line). 



nitride of this process can be only poorly constrained due 
to the huge spread in the observed metallicities of GCs. 
Owing to the variety of processes contributing to the mix- 
ing of metals with the shell, the spread in metallicity can 
reflect, at least in part, a cluster-to-cluster variation of 7. 

It is also worth noting that the evolution of the metal- 
licity [Fe/H] as a function of My is almost flat. This is in 
agreement with the fact that no evident trend is shown in 
the observed My vs. [Fe/H] relationship. 

A similar comparison of the results of the models with 
the observations can be made taking into consideration 
the masses of GCs instead of the magnitudes. In this case, 
we do not have to use eq. (31) but we can directly com- 
pare the swept-up mass at the moment of the onset of 
the gravitational instabilities with the mass tabulated by 
Pryor & Meylan (1993) (see Fig. 4). The results are plot- 
ted in Fig. 10. The conclusions we can draw from this plot 
are the same: a very low mixing efficiency can be ruled 
out and, in order to reproduce the GCs with the largest 
metallicities, we have to assume that more than 50 % of 
the metals mix with the surrounding shell in a time-scale 
shorter than the time needed for the onset of gravitational 
instabilities. 

4.2.2. The Magnitude-Metallicity relation as a function 

of mi 

Many authors believe that the IMF of primordial objects 
should have been top-heavy, owing to the fact that the 
CMB radiation keeps the gas warm enough to raise the 
Jeans mass. We therefore varied the lower limit of the 
IMF distribution to; in order to see what effect an IMF 




Fig. 11. Same as Fig. 9 but with varying lower mass 
to;. This parameter ranges the values 0.1 (solid line), 0.5 
(dotted line), 1 (dashed line), 3 (long-dashed line) and 10 
(dot-dashed line). 

biased towards massive stars can have in the metallicity 
of GCs. 

The Magnitude-Metallicity relation as a function of m; 
is shown in Fig. 11. In this set of models, the thermaliza- 
tion efficiency 77 and the mixing efficiency 7 are set to be 
10 -4 and 0.5, respectively. As we can see in this figure, an 
IMF strongly biased towards massive stars (to; = 10) can 
justify the most metal-rich GC observed. If the IMF is so 
strongly top-heavy, mixing efficiencies larger that 0.5 can- 
not be allowed. This mixing efficiency is consistent with 
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Fig. 12. Same as Fig. 9 but with varying thermalization 
efficiency 77. This parameter ranges the values 5 • 1CP 3 
(solid line), 1CT 3 (dotted line), 10~ 4 (dashed line), 1(T 5 
(long-dashed line) and 10 -6 (dot-dashed line). 



Fig. 13. Same as Fig. 9 but with varying external pres- 
sure P4. This parameter ranges the values 1 (solid line), 5 
(dotted line), 10 (dashed line), 50 (long-dashed line) and 
100 (dot-dashed line). 



the value found by Recchi et al. (2001) in their simulations 
of IZwl8. Owing to the fact that the selected sample of 
GCs is not perfectly coeval, small cluster-to-cluster vari- 
ations of mi are possible, explaining therefore, at least in 
part, the observed scatter. 

Also for these models, [Fe/H] does not vary so much 
with My. Only for models with large mi a correlation 
arises, in the sense that the faintest (and less massive) 
GCs are the most metal-rich. This is due to the fact that 
t* grows with M§. At large t, Mz is almost independent 
on t (see Fig. 6), whereas M sw grows linearly with t (eq. 
21). From eq. (27) we can see that Z should decrease. This 
trend has been found also by P99. Statistical properties 
of a selected sample of old GCs (the so-called Old Halo 
GCs) seem to confirm this trend (Parmentier & Gilmore 
2001), although the scatter is huge. 

4.2.3. The Magnitude-Metallicity relation as a function 
of r\ 

The thermalization efficiency r\ is the ratio of the mechan- 
ical energy of a single SN to the internal energy gained 
by the ISM after the explosion. As we have seen in Sect. 
3.2, the value of 77 can be estimated analytically (eq. 9) 
but, due to the large uncertainties in this evaluation, we 
decided to keep the thermalization efficiency as a free pa- 
rameter. However, owing to the very large densities in the 
core of PGCs, this value should be very low. We therefore 
decided to span the parameter range [5 • 10~ 3 , 10~ 6 ]. 

The Magnitude-Metallicity relation as a function of 77 
is shown in Fig. 12. In this set of models, the lower limit 
of the IMF distribution mi and the mixing efficiency 7 
are set to be 1 and 0.5, respectively. As we can see in this 
figure, below 77 = 10~ 4 , the dependence of the metallicity 



with 77 is very weak. Above this value, the dependence 
becomes stronger due to the non-linear dependence of t* 
(and consequently Mz(t*) and M sw (t*)) on 77 (see eq. 26). 

4.2.4. The Magnitude-Metallicity relation as a function 
of P 4 

So far we have analyzed the dependence of the My vs. 
[Fc/H] relation as a function of internal parameters of the 
PGC. Of course also external parameters play a role. In 
this subsection we analyze the evolution of the metallicity 
of GCs as a function of the external pressure P4. This 
dependence is shown in Fig. 13. 

The effect of the external pressure is similar to the ef- 
fect of the thermalization efficiency 77, in the sense that 
above P4 — 10 this effect is negligible and the tracks lie 
close to each other. For low values of the external pres- 
sure, strongly non-linear effects on the eq. (26) create the 
knee observed for the model with P4 = 1 (solid line) and 
P4 = 5 (dotted line). Smaller values of P4 violate the con- 
dition 3c s Vb < V2Gnporl (see Sect. 3.5) and cannot be 
taken into account in this simplified approach. As stated 
in Sect. 4.1, these low pressure environments are unlikely 
in protogalactic backgrounds. 

5. Conclusions 

We performed an analytical study of the evolution of a 
proto-globular cluster (PGC) assuming an initial burst 
of star formation occurring inside the core radius of the 
initial gaseous distribution. We followed the evolution of 
the shock wave formed after energy release from the first 
SNell and we considered the possibility that gravitational 
instabilities may arise in the swept-up shell, leading to the 
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formation of a second population of stars, with metallicity 
larger than zero, owing to the pollution from the ejecta of 
the first supernovae. Our results can be summarized as 
follows: 

— Is it possible to calculate a lower mass of the PGC 
assuming that the onset of the gravitational instability 
able to produce the second generation of stars should 
occur before the shell reaches the edge of the PGC 
gaseous distribution. With reasonable choices of the 
parameters, this threshold mass is of the order of 10 6 
M Q , in agreement with other similar studies (SWOO, 
Brown et al. 1995). 

— In order to produce a second stellar population with a 
metallicity in agreement with the one observed in the 
most metal-poor GCs, it is necessary to mix at least 
5 - 10 % of the metals produced in the central burst 
of star formation with the surrounding swcpt-up shell. 
The most metal-rich GCs are consistent with a very 
large mixing efficiency (larger than 0.5; i.e. more than 
50% of the metals mix with the surrounding shell). 

— The evolution of theoretical tracks in the My vs. 
[Fe/H] plot is almost constant. This is consistent with 
the fact that no obvious trend can be found in the ob- 
served My vs. [Fe/H] relationship. The scatter of the 
metallicity is therefore due mostly to scatter in the 
internal parameters of the cluster. 

— A top-heavy IMF (with lower mass mi larger than 1 
M Q ) seems to be required in order to explain the GCs 
with the largest mctallicities. 

— The huge spread in the metallicity observed in the sam- 
ple of GCs can be explained, at least in part, with 
cluster-to-cluster variations of the structural parame- 
ters 7, r\ and m;, as well as variations of the external 
pressure P4. 
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